DETERMINING PROJECTIONS AND FUNCTIONALS FOR 
WEAK SOLUTIONS OF THE NAVIER-STOKES EQUATIONS 



MICHAEL HOLST AND EDRISS TITI 

Abstract. In this paper we prove that an operator which projects weak solutions of 
the two- or three-dimensional Navier-Stokes equations onto a finite-dimensional space 
is determining if it annihilates the difference of two "nearby" weak solutions asymptot- 
ically, and if it satisfies a single appoximation inequality. We then apply this result to 
show that the long-time behavior of weak solutions to the Navier-Stokes equations, in 
both two- and three-dimensions, is determined by the long-time behavior of a finite set 
of bounded linear functionals. These functionals are constructed by local surface aver- 
ages of solutions over certain simplex volume elements, and are therefore well-defined 
for weak solutions. Moreover, these functionals define a projection operator which satis- 
fies the necessary approximation inequality for our theory. We use the general theory to 
establish lower bounds on the simplex diameters in both two- and three-dimensions. Fur- 
thermore, in the three dimensional case we make a connection between their diameters 
and the Kolmogoroff dissipation small scale in turbulent flows. 
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1. Introduction 

Consider a viscous incompressible fluid in £1 C IR d , where Q is an open bounded 
domain with Lipschitz continuous boundary, and where d — 2 or d — 3. Given the 
kinematic viscosity v > 0, and the vector volume force function f(x,t) for each x E £1 
and t e (0, oo), the governing Navier-Stokes equations for the fluid velocity vector u = 
u(x, t) and the scalar pressure field p = p(x, t) are: 

Ou 

— - uAu + ■ V)u + Vp = f inQx(0, oo), (1.1) 

\J v 

V-m = in fix (0,oo). (1.2) 
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Also provided are initial conditions u(Q) = u , as well as appropriate boundary condi- 
tions on <9f2 x (0, oo). 

The notion of determining modes for the Navier-Stokes equations was first introduced 
in [13] as an attempt to identify and estimate the number of degrees of freedom in turbu- 
lent flows (cf. [9] for a thorough discussion of the role of determining sets in turbulence 
theory). This concept later led to the notion of Inertial Manifolds [14]. An estimate of the 
number of determining modes was given in [12] and later improved in [21]. The notion 
of determining nodes, and other more general determining concepts, were introduced 
in [15]. In [16] the notion of determining nodes was discussed in detail, and estimates 
for their number were reported in [20], and later improved in [21]. In [17] (see also [19]) 
the concept of determining volume elements was presented, and a connection was estab- 
lished between this concept and Inertial Manifolds. A generalized and unified theory of 
all of the above was recently presented in [5, 6]. 

Bounds on the number of determining modes, nodes, and volumes are usually phrased 
in terms of a generalized Grashof number, which is defined for the two-dimensional 
Navier-Stokes equations as: 



where Ai is the smallest eigenvalue of the Stokes operator and p = yAi is the related 
(best) Poincare constant. Here, F = limsup t ^ 00 (/ Q \f(x, i)] 2 ) 1 / 2 if / G L 2 (Q) for 
almost every t, or F = limsup^^ \Ai||/||i?-i(n) if / € H^ 1 ^) for almost every t. 

The best known estimate for the determining set size for the two-dimensional Navier- 
Stokes equations with periodic boundary conditions and i/ 2 -regular solutions is of order 
Gr [21]. In obtaining their estimate, the authors relied on the fact that the domain had no 
physical boundaries to shed vorticity, which made available some convenient properties 
of H 2 -regular solutions. However, in the two-dimensional case with no- slip boundary 
conditions, to our knowledge the best estimate on the cardinal of any determining set 
(modes, nodes, or volumes) that can be obtained is of order Gr 2 , even for if 2 -regular 
solutions. 

Due to the Sobolev Imbedding Theorem H 2 ^ C° (which holds in dimensions 1, 

2, and 3), or rather due to the failure of the imbedding H 1 C° in dimensions 2 and 

3, determining node analysis is necessarily restricted to if 2 -regular solutions to make 
sense of point- wise values. However, when talking about determining modes or volume 
elements, it is sufficient for functions to be H 1 -regular, so that these concepts also make 
sense for weaker solutions. To construct a general analysis framework for the case of 
weak H 1 solutions, we can begin by defining notions of determining projections and 
determining junctionals for weak solutions. (The standard spaces H, V, and V are 
reviewed fully in §2.) 

Definition 1.1. Let f(t),g(t) 6 V be any two forcing functions satisfying 

lim \\f(t)-g(t)\\ v , = 0, (1.3) 

t— >oo 

and let u,v E V be corresponding weak solutions to (1.1)— (1.2). The projection operator 
Rn '■ V (-> Vn C L 2 {VL), N = dim(Vjv) < oo, is called a determining projection for 
weak solutions of the d-dimensional Navier-Stokes equations if 

lim \\R N (u(t)-v(t))\\ L 2 {Q) = 0, (1.4) 

implies that 

lim \\u(t) -v(t)\\ H = 0. (1.5) 

t— >oo 
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Given a basis {(j)i}f =1 for the finite-dimensional space V N , and a set of bounded linear 
functionals {k}^ from V, we can construct a projection operator as: 

N 

RnU = y^Jjju)^. (1.6) 

i=l 

The assumption (1.4) is then implied by: 

lim \k(u(t) -v(t))\ = 0, i = l,...,iV (1.7) 

t— >oo 

so that we can ask equivalently whether the set {Z i }^ 1 forms a set of determining func- 
tionals (see [5, 6]). The analysis of whether or {£j}|Li are determining can be reduced 
to an analysis of the approximation properties of R^. Note that in this construction, the 
basis {4>i]f = i need not span a subspace of the solution space V, so that the functions 0; 
need not be divergence-free for example. Note that Definition 1.1 encompasses each of 
the notions of determining modes, nodes, and volumes by making particular choices for 
the sets of functions {(j>i}f =1 and {k}f =l (see [19, 20]). 

In this paper, we will employ Definition 1.1 to extend the results of [5, 6] to the more 
general setting of ii/ -1 -regular solutions. In particular, we will show that if a projection 
operator R N : V i-> Vjv C L 2 (Vi), N = dim(Vjv) < oo, satisfies an approximation 
inequality for 7 > of the form, 

||u - R n u\\l^(q) < CiN^\\u\\ H i {n) , (1.8) 

then the operator Rn is a determining projection in the sense of Definition 1.1, provided 
N is large enough. We will also derive explicit bounds on the dimension N which guar- 
antees that R N is determining. While we gain generality in our approach here, we also 
lose something in the balance: the bounds obtained here are generally of order Gr 2 , 
whereas the bounds in [5, 6] (requiring H 2 -regularity) are of order Gr. 

Outline of the paper. Preliminary material is presented in §2, including some inequali- 
ties for bounding the nonlinear term appearing in weak formulations of the Navier-Stokes 
equations. In §3, a finite element interpolant due to Scott and Zhang is presented, which 
(unlike nodal interpolation) is well-defined for H 1 -functions. It is shown that the inter- 
polant satisfies the approximation assumption (1.8) for f/^-functions on arbitrary poly- 
hedral domains in both two and three dimensions; most of the details are relegated to the 
Appendix. In §4, we consider the two-dimensional Navier-Stokes equations, and derive 
bounds on the dimension iV of the space Vn, employing only the approximation assump- 
tion (1.8). As an application of this general result, we employ some standard assumptions 
about simplex triangulations of the domain (discussed in §3) and derive lower bounds on 
the simplex diameters, sufficient to ensure that the SZ-interpolant is a determining projec- 
tion (equivalently, that the simplex surface integrals forming SZ-interpolant coefficients 
are a determining set of linear functionals). We extend these results to three dimensions 
in §5, by requiring (following [7]) that weak solutions satisfy an additional technical as- 
sumption (due to the lack of appropriate global a priori estimates), which is related to 
the natural notion of mean dissipation rate of energy. 

2. PRELIMINARY MATERIAL 

We briefly review some background material following the notation of [8, 23, 25, 26]. 
Let Cl C IR d denote an open bounded set. The imbedding results we will need are 
known to hold for example if the domain Cl has a locally Lipschitz boundary, denoted 
as fi G C 0,1 (cf. [1]). For example, open bounded convex sets flcK 1 * satisfy f2 e C ' 1 
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(Corollary 1.2.2.3 in [18]), so that convex polyhedral domains (which we consider here) 
are in C 0,1 . 

Let H k (Q) denote the usual Sobolev spaces W k,2 {Vt). Employing multi-index nota- 
tion, the distributional partial derivative of order |a| is denoted D a , so that the (integer- 
order) norms and semi-norms in H k (n) may be denoted 

k 

Nljr*(n) = Yl M^Mikn)' \ u \ 2 HHn) = Yl IPNUw < j < /c, 

3=0 \a\=j 

where \ represents the measure of fi. Fractional order Sobolev spaces and norms may 
be defined for example through Fourier transform and extension theorems, or through 
interpolation. A fundamentally important subspace is the k = 1 case of 

H k (£l) = closure of in H k (tt), 

in which the Poincare Inequality reduces to: If is bounded, then 

IMIi^n) < p{fy\u\m(n), VueH%(Q). (2.1) 

The spaces above extend naturally (cf. [25]) to product spaces of vector functions u = 
(u\, u 2 , ■ ■ ■ , Ud), which are denoted with the same letters but in bold-face; for example, 
Hp(fi) = (//^(fi)) . The inner-products and norms in these product spaces are extended 
in the natural Euclidean way; the convention here will be to subscript these extended 
vector norms the same as the scalar case. 

Define now the space V of divergence free C°° vector functions with compact support 

as 

V = {0G Cg°(ft) | V-0 = O}. 

The following two subspaces of L 2 (f2) and Hj(f2) are fundamental to the study of the 
Navier-Stokes equations. 

H = closure of V in L 2 (f2), V = closure of V in Hq(^). 

To simplify the notation, it is common (cf. [8, 25]) to use the following notation for 
inner-products and norms in H and V: 

(u,v) = {u,v) H , \u\ = \\u\\ H , {{u,v)) = {u,v) v , \\u\\ = \\u\\ v . (2.2) 

The Navier-stokes equations (1 .1)— (1.2) are equivalent to the functional differential 
equation: 

du 

— + uAu + B(u,u) = f, u(0)=u . (2.3) 
at 

The Stokes operator A and bilinear form B are defined as 

Au = -PAu, B(u, v) = P[(u ■ V>], 

where the operator P is the Leray orthogonal projector, P : Hq i— > V and P : L 2 i-> H, 
respectively. 

Weak formulations, which we consider shortly, will use the bilinear Dirichlet form 
((•, •)) and trilinear form &(•, •, •) as: 

((■u, v)) = (Vw, Vi>), b(u, v, w) = (B(u, v),w) = (P((u ■ V)v),w). 

(Note that thanks to the Poincare inequality (2.1), the form ((•,•)) ^ s actually an inner- 
product on V, and the induced norm || • || = ((•, -)) 1 / 2 is in fact a norm on V, equivalent 
to the iJ 1 -norm.) 
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A priori bounds can be derived for the form &(•, ■) (cf. [8, 22, 25]). In particular, if 
Vl C W 1 , then the trilinear form b(u, v, w) is bounded on V x V x V as follows: 

d = 2 : \b(u,v,w)\ < 1| || | i/^^) 1^ I j^t^ || ^ || J^^^^ | ^ | i/? c ^) , (2.4) 

d = 3 : \b{u,v,w)\ < 2||«||^ n) |tt|^ 4 (n) |«| ff i(n)||w||^ 4 (n) klfff(n)- ( 2 - 5 ) 
Moreover, from Holder inequalities we have for d = 2 or d = 3: 

\b(v,u,v)\ < ||VM|U<»(n)||u||za(n)- ( 2 - 6 ) 

3. Polynomial interpolation in Hj(fi) 

An example of a projection operator which satisfies the approximation assumption (1.8) 
is that used for defining determining volumes [19]; we examine now powerful alternative 
operator. Let Q C IR d be a d-dimensional polygon, exactly triangulated by (for example) 
Delaunay triangulation [11], with quasi-uniform, shape-regular simplices, the vertices of 
which will form a set of iV generalized interpolation points in our analysis. Note that 
for quasi-uniform, shape-regular triangulations in W 1 (see [4] for detailed discussions), it 
holds that 

C \Q\h- d < N <C' Q \Q\h~ d , (3.1) 

where h is the maximum of the diameters of the simplices, and where Co and C' Q are 
universal constants, independent of both N and h. The parameter h will be referred to 
as the characteristic parameter, or characteristic length scale, of such a quasi-uniform 
shape-regular mesh. 

It should be noted that given some initial triangulation satisfying (3.1), repeated bi- 
section [2] or octa-section [27] (quadra-section in 2D) of each simplex can be performed 
in such a way as to guarantee non-degeneracy asymptotically, in that the quasi-uniformity 
and shape-regularity are preserved. Therefore, inequality (3.1) can be made to hold, for 
the same universal constants, for finer and finer meshes in a nested sequence of simplex 
triangulations. 

To properly define a continuous piecewise-linear nodal interpolant of a function u G 
H x {yi) based on the nodes of a triangulation of Vt, the particular function u must be 
bounded point-wise. This will be true if the function u is continuous in Vt, hence uni- 
formly continuous on Vt. One of the Sobolev imbedding results (cf. [1]) states that if 
Vt C W 1 satisfies £1 G C 0,1 , then for nonnegative real numbers k and s it holds that 
H k (fl) <— > C s (Cl), k > s + |. This implies that for d — 1, the interpolant can be 
correctly defined, since H 1 ^) is continuously imbedded in C°(f2). However, in higher 
dimensions, H 1+a (Q) ^ C°(Q) only if a > when d = 2, or if a > 1/2 when d = 3. 
While it may be possible to use the nodal interpolant and a regularity assumption such as 
u E H 1+a (Q) for appropriate a > 0, an alternative approach is taken here. 

The generalized interpolant due to Scott and Zhang [24] can be defined for if 1 -functions 
in both two and three spatial dimensions. The SZ-interpolant I h is constructed from a 
combination linear interpolation and local averaging on faces and edges of simplices, 
and has optimal approximation properties even in the case of H 1 -functions. 

Lemma 3.1. For the SZ-interpolant of u G H 1 ) +a (£l), a > 0, it holds that 

\\u - I h u\\tf(n) < Cih 1+a \u\ H i+ a{n) . 



Proof. See the appendix for a condensed proof following [3, 24]. 



□ 
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Note that both the usual nodal interpolant and the SZ-interpolant I h can be written as 
a linear combination of linear functionals: 

N 

I h u(x) = s ^2<p i (x)l i (u). 

1=1 

In either case, the set of functions {<pi}f =l is the usual continuous piecewise-polynomial 
nodal finite element basis defined over the simplicial mesh, satisfying the Lagrange prop- 
erty at the vertices of the mesh: 

4>i(xj) = 8ij. 

The difference between the two interpolants is simply the choice of the linear functionals: 
in the case of the nodal interpolant, the functionals are delta functions centered at the 
vertices of the mesh; in the case of the SZ-interpolant, they are defined in terms of a 
bi-orthogonal dual basis (see the Appendix). 



4. The Two-dimensional Navier-Stokes Equations 

A general weak formulation of the Navier-Stokes equations (1.1)— (1.2) can be written 
as (cf. [8, 25]): 

Definition 4.1. Given / G L 2 ([0, T]; V), a weak solution of the Navier-Stokes equations 
satisfies u G L 2 ([0,T); V) n C w ([0,T);H), du/dt G Ll oc ((0,T);V'), and 

du 

< — , v > +v((u, v)) + b(u, u, v) =< f, v >, Vf G V, for almost every t, (4.1) 

u(0) = u . (4.2) 

Here, the space C w ([0, T\ \ H) is the subspace of L°°([0, T]; H) of weakly continuous 
functions, and < •, • > denotes the duality pairing between V and V, where H is the 
Riesz-identified pivot space in the Gelfand triple V C H = H' C V '. Note that since 
the Stokes operator can be uniquely extended to A : V h->- V' , and since it can be shown 
that B : V x V M- V (cf. [8, 26] for both results), the functional form (2.3) still makes 
sense for weak solutions, and the total operator represents a mapping V H V. 

In the two-dimensional case, for a forcing function / G L°°([0,T]; V), there exists 
a unique weak solution u G L 2 ([0,T]; V) n C w ([0,T);H) (cf. [8, 26]). Consider now 
two forcing functions f,g G L 2 ([0, oo]; V) and corresponding weak solutions u and v 
to (2.3) in either the two- or three-dimensional case. Subtracting the equations (2.3) for 
u and v yields an equation for the difference function w = u — v, namely 

^ + vAw + B(u,u) - B(v,v) = f - g. (4.3) 
dt 

Since the residual of equation (4.3) lies in the dual space V, for almost every t, we can 
consider the dual pairing of each side (4.3) with a function in V, and in particular with 
w G V, which yields 

< ~^jr^ w > +HIHI 2 + H u i u i w ) ~ H v i v i w ) =< f ~ 9i w > f or almost every t. 
It can be shown (cf. [25], Chapter 3, Lemma 1.2) that 
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in the distribution sense. It can also be shown [8, 25] that b(u,v,w) = —b(u,w,v), 
Vu,v,w G V, so that b(w,u,w) = b(u,u,w) — b(v,v,w). Therefore, the function 
w = u — v must satisfy 

-^|w| 2 + z/||w|| 2 + b(w, u, w) =< f — g,w > . (4.4) 

The following generalized Gronwall inequality will be a key tool in the analysis to 
follow (see [12] and [19]). 

Lemma 4.2. Let T > be fixed, and let a(t) and /3(t) be locally integrable and real- 
valued on (0, oo), satisfying: 

Y pt+T j rt+T 



liminf— / a{r)dr = m > 0, lim sup — / a (r)dr = M < oo, 



i r t+T 

lim- /3+(r)dr = 0, 

t— >oo 1 J. 



where a = max{— a,0} and (3 + = max{/3,0}. If y(t) is an absolutely continuous 
non-negative function on (0, oo), and y(t) satisfies the following differential inequality: 

y'(t) + a{t)y{t) < (3{t), a.e. on (0, oo), 

then lim^oo y(t) = 0. 

The main two-dimensional results are now given; we assume that f2 C IR 2 is an open 
bounded domain with Lipschitz continuous boundary. 

Theorem 4.3. Let f(t),g(t) G V be any two forcing functions satisfying 

}im\\f(t)-g(t)\\ v , = 0, 

and let u,v G V be the corresponding weak solutions to (1. !)-(!. 2) for d = 2. If there 
exists a projection operator Rn : 7 4 Vjv C L 2 (Q), N = dim(VAr), satisfying 

lim \\R N (u(t) - v(t))\\ L 2 (n) = 0, 
and satisfying for 7 > the approximation inequality 

||tt - RnuWl^q) < CW _7 IMI//i(n), 

then 

lim \u(t) -v(t)\ = 

t— >OD 

holds if N is such that 

1 

1 . \^ 



00 > N > C ( — lim sup 



V 



where C is a constant independent of v and f. 

Proof. Using the notation (2.2), we begin with equation (4.4), employing the inequal- 
ity (2.4) along with Cauchy-Schwarz and Young's inequalities to yield 

\w\ 2 + ^||u>|| 2 < \\u\\ \w\ \\w\\ + 11/ — 



2dV 



< -\\uf\w\ 2 + -\\f - gf v , + ^\\w\ 
v v 2 
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Equivalently, this is 

^-\w\ 2 + u\\w\\ 2 - -|M| 2 M 2 < -11/ - g\\v>- 

at v v 

To bound the second term on the left from below, we employ the approximation assump- 
tion on R N , or rather the following inequality which follows from it: 

|2 



\W\ 



<2iV- 2 ^ 1 2 || W || 2 + 2||i? i v U ;||i 2(f , ) , 



which yields 



d , l2 fuN 2 ^ 2„ ll2 \ , l2 2.., Il2 vN 2 ~* „ „ ll2 

M + -^r - -Ml M < -11/ - ^ v' + -^HI^HI 2 



This is of the form 

— to + a w < p, 
at 

with obvious definition of a and p\ 

The generalized Gronwall Lemma 4.2 can now be applied. Recall that both ||/ — 
g\\v — > and ||-RjvHU 2 (!r2) — > as t — > oo by assumption. Since it is assumed that 
u and v, and hence w, are in V, so that all other terms appearing in a and (3 remain 
bounded, it must hold that 

i r t+T i r t+T 

lim — / /3 + (r)dr = 0, limsup— / a (r)dr < oo. 
It remains to verify that for some fixed T > 0, 

-t+T 



1 f t+1 

limsup — / a{r)dr > 0. 



This means we must verify the following inequality for some fixed T > 0: 

iv27 > m c limsup i r +T ?jh^ \ = m limsup i /" t+T ii^p^ (4 5) 

V \ t^oo T J t V ) V 1 t-^oo T J t 

The following a priori bound on any weak solution can be shown to hold (this is a simple 
generalization to / € V of the bound in [8] for / E H): 

-t+T 



1 f + 2 
limsup- / \\u(T)\\ 2 dT < —limsup \\f(t)\\y,, 

t— >oo J- Jt v t-^oo 



for T = p 2 jv > 0, where p is the best constant from the Poincare inequality (2.1). 
Therefore, if 

N 2y > 8C 2 ("J_li msU p|| /( t)|| y ,y > M (Alimsup||/(t)||^) , (4.6) 

implying that (4.5) holds, then by the Gronwall Lemma 4.2, it follows that 

lim \w(t)\ = lim \u(t) - v(t)\ = 0. 

t— >oo t— >oo 

□ 

Assume now that Q C IR 2 is also polyhedral, and can be exactly triangulated with 
a quasi-uniform, shape-regular set of simplices of maximal diameter h = 0(N' 1 ^ 2 ), 
where iV is the number of vertices in the triangulation (see §3). As an application of the 
general result above, we establish a lower bound on the simplex diameters of such a trian- 
gulation, which ensures that the SZ-interpolant is a determining projection (equivalently, 
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that the simplex surface integrals forming SZ-interpolant coefficients are a determining 
set of linear functionals). 

Corollary 4.4. The SZ-interpolant is determining for the two-dimensional Navier-Stokes 
equations if the diameter h of the simplices is small enough so that 

oo > bT 2 > C { — limsup 

Proof. Since h = 0(N~ l l 2 ) for quasi-uniform, shape-regular triangulations in two di- 
mensions, taking a = in Lemma 3.1 yields 

\\u - I h u\\jj»(n) < Cih\u\ H i(u) < CiN~ 1/2 \\u\\ H i {n) . 

Therefore, the SZ-interpolant I h satisfies the approximation inequality (1.8) for 7 = 1/2. 
The corollary then follows by application of Theorem 4.3. □ 

Remark 4.5. If f G H, then we have in fact a strong solution, i.e. u G H 2 (Q), and 
the interpolation Lemma 3.1 may be applied with a = 1. This falls into the theoretical 
framework of [5, 6], and in the periodic case they have shown that N « Gr, whereas the 
above result for the no-slip case states that N w Gr 2 . Whether the no-slip case may be 
improved to N « Gr with additional regularity (f G H) is unclear, due to the lack of an 
analogous identity to 

(B(w,w),Aw) = 0, 

which holds for the two-dimensional periodic case. In physical terms, in two dimen- 
sions this identity illustrates the lack of a boundary vorticity shedding source when the 
boundary is absent. 

5. The Three-dimensional Navier-Stokes Equations 

The lack of appropriate a priori estimates in the three-dimensional case requires a 
modification of the approach taken for the two-dimensional case in the previous section. 
However, the interpolation results we have employed are dimension-independent, and 
by following the analysis approach of [7] very closely, we can obtain similar results for 
the three-dimensional case. Again we require only that / G V, but we also assume the 
existence of a unique weak solution to the three-dimensional Navier-Stokes equations. 
An additional technical assumption is that some measure of the mean rate of energy 
dissipation be finite, namely: 

v f t+T u , 
= inf lim sup — / Vw L, dr < oo. 
t>o T J t 

This assumption implies that eventually the weak solution for the three-dimensional 
Navier-Stokes equations becomes unique, and also in the case / G H the weak solu- 
tion eventually becomes strong. But this assumption does not imply anything about the 
transients, since the quantity is required to be finite only for large time. We assume again 
that n C M 3 is an open bounded domain with Lipschitz continuous boundary. 

Theorem 5.1. Let f(t),g(t) G V be any two forcing functions satisfying 

lim ||/(*)-<7(t)||y = 0, 

t— >oo 

and let u,v G V be the corresponding weak solutions to (!.!)-( 1.2) for d = 3. If there 
exists a projection operator R N : V t-> Vat C L 2 (Q), N = dim(VTv), satisfying 

lim \\R N (u(t) - v(t))\\ L2{u) = 0, 




10 M. HOLST AND E. S. TITI 

and satisfying for 7 > the approximation inequality 

\\u - Rnu\\ L 2( Q) < CiiV~ 7 ||u||Hi(n), 

then 

lim \u(t) -v(t)\ = 

t— >oo 

holds if N is such that 



00 > N > C (- inf (limsup ^ I \\Vu(s)\\ Laam dsX\ 

\U T>0 [ 1 J f J J 

where C is a constant independent of v, f, and u. 

Proof. Beginning with equation (4.4), the inequality (2.6) is employed along with Cauchy- 
Schwarz and Young's inequalities to yield 

^ \w\ 2 + ^||u>|| 2 < || Vw||l° c (^)|u>| 2 + ||/ — gr||y/||iy|| 



2dt 



< ||Vm|| l ^M 2 + —11/ - g\\ 2 v , + ^|H| 2 



Equivalently, 

^t\w\ 2 + is\\w\\ 2 - ||Vm||l-»(c)|w| 2 < -||/ - g\\v>- 

dt v 
To bound the second term on the left from below, we employ a consequence of the 
approximation assumption on R]y, namely the inequality 



\w 



2 <2N- 2 ^C 2 1 \\w\\ 2 + 2\\R N w\\ 2 L2{n) , 



which yields 



M + ^2" " V« L « hi" < "11/ " ^ll 2 /' + ^2-p7V^|| 2 



dt 1 1 V^Cf" ~ l|L / 11 u lu ~ llv ' ~(^» £Lfru, »i< , (ny 
This has the form 

— u> + a kw < p, 
at 

with again obvious definition of a and (3. 

The analysis now proceeds exactly as in the proof of Theorem 4.3, so that all that 
remains is to check again that for some fixed T > 0, 

-t+T 



1 f t+1 

limsup — / a(r)dr > 0. 

t— >oo J- Jt 



Thus, we must prove our assumption on N guarantees for a fixed T > that 

2C 2 1 f t+T 

N 2j > — - lim sup — / \\Vu\\ L oo {n) dT. (5.1) 

V t-*oc J- Jt 

If we select T* > such that 

/ 1 rt+T x j ,i+T* 

2mf limsup — / || Vu(s)\\ L ^^ Q) ds ) > limsup — / \\Vu(s)\\ L oo^ds, 

T >° V t->oo 1 Jt J t^oo 1* Jt 

then our assumption gives 

AC 2 ( 1 f t+T * \ 

N ^ > _L in f ii msU p / ||Vu(s)|| L ~ (n )ds (5.2) 

V T»>0 y f^oo 1* J t J 

which implies (5.1). The theorem then follows by the Gronwall Lemma 4.2. □ 
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Assume now that Q C IR 3 is also polyhedral, and can be exactly triangulated with 
a quasi-uniform, shape-regular set of simplices of maximal diameter h = 0(A r_1 / 3 ), 
where N is the number of vertices in the triangulation. As an application of the general 
three-dimensional result above, we will establish a lower bound on the simplex diameters 
of such a triangulation, which ensures that the SZ-interpolant is a determining projection 
(and that the simplex surface integrals forming SZ-interpolant coefficients are a deter- 
mining set of linear functionals). 

Corollary 5.2. The SZ-interpolant is determining for the three-dimensional Navier-Stokes 
equations if the diameter h of the simplices is small enough so that 

oo > h~ 2 > C ( - inf < limsup — f \\^u(s)\\Loo^ds\ J . 

Proof Since h = O^N- 1 / 3 ) for quasi-uniform, shape-regular triangulations in three 
dimensions, taking a = in Lemma 3.1 yields 

\\u - i h u\\ L 2 (n) < Ci/i|«|#i( n ) < CiN^ 1/3 \\u\\ H i {n) . 

Therefore, the SZ-interpolant Ih satisfies the approximation inequality (1.8) for 7 = 1/3. 
The corollary then follows by application of Theorem 5.1. □ 

Appendix: Approximability of the Scott-Zhang Interpolant 

We will sketch the proof of the approximability result for the SZ-interpolant given as 
Lemma 3.1; we will follow quite closely the proof given in [3, 24]. As throughout this 
paper, we assume that G C 01 , and that the given exact simplicial triangulation of Vl is 
both shape-regular and quasi-uniform. 

The proof of Lemma 3 . 1 will follow easily from the following result (see the comments 
at the end of this appendix). 

Lemma 5.3. For the SZ-interpolant ofu G Hl +a {VL), a > 0, it holds that 

\\u - I h u\\ L 2 {n) < Ci^ 1+a M//i+«(n)- 

To prove Lemma 5.3, we will begin by defining carefully the SZ-interpolant. Let Th = 
{ T i}f=i be the given quasi-uniform, shape-regular mesh of d-simplices which exactly 
triangulate the underlying domain fi, and let Vl h = {xi]f =1 be the set of vertices of these 
d- simplices. Define 

V h = span{<j )l (x)}l 1 cH 1 (n), 
where {0j(x)} is the set of standard continuous piecewise linear (nodal) basis functions. 
The nodal basis satisfies the Lagrange relationship at the vertices (which are exactly the 
"nodes" in this setting): 

4>i{xj) = S i:j . 

Now, for each vertex we select (arbitrarily) an associated (d — l)-simplex cr, from the 
given simplicial mesh satisfying only: 

(1) Xi G <7i, and (2) ai C Oil if Xi G <9f2. 
In other words, for a given vertex x« we pick an arbitrary (d — l)-simplex from edges or 
faces of the (i-simplices which contain Xj as a vertex. In two-dimensions, we are picking 
the edge of one of the triangles that have Xi as a vertex; in three-dimensions, we are 
picking the face of one of the tetrahedra which have Xi as a vertex. The only restriction 
on this choice is near the boundary: if x^ is on the boundary, then the [d — l)-simplex 
we pick must be one of the edges or faces of the a simplex which lies exactly on the 
boundary (such a choice is always possible). 
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In each (d — l)-simplex <7i, we number the generating vertex Xj first in the set of 
vertices of a h denoted {x iy j}^ =1 . (I.e., we set x iy i = Xj.) For each <jj, we also have a 
(d — 1) -dimensional nodal basis {(f>i,j}j=i, where again we set 0^1 = 0j. There exists an 
associated L 2 (crj)-dual (bi-orthogonal) basis {if)ij} satisfying 

/ tfj id (x)(f)i : k(x)dx = 5jk, j,k = l,...,d. 

Again we take ip it i = ipi, Vxj G f4 Note that ^ and <pj also satisfy a bi-orthogonal 
relationship, namely J ipi4>jdx = 0, % ^ j. We define now the SZ-interpolant as 

N 

I h : H\n) ^ ^(O), J fcU (ar) = ^^(x)^(w), Z^u) = / V;(£M£)cf£. 

i=i * 

Thanks to the Trace Theorem [1], the interpolant Ihu(x) is well-defined at nodal values 
even for u G H l {Vt), since >■ L?{o.i). Almost by construction, one can show [24] 

that 

• 4 '■ H 1 ^) i — ^ Vft(Q) is a projection 

where Vo^ is the subset of having zero trace on the boundary of O. Thus, 4 P re ~ 
serves homogeneous Dirichlet boundary conditions. Using homogeneity arguments, the 
following stability result for the interpolant is established in [24] . 

Lemma 5.4. For any r G 4. if the support region of r is defined as the set S T = 

interior (U{fj | D f =fi 0, 7$ G 4}), ^en if fa?k£s 

||4w|k-(r) < C^/i fe - m |M| Hfc(Sr) , 0<m<Z, f > 1/2. 

fc=0 

Proof. See the proof of Theorem 3.1 in [24]. □ 

The proof of the Scott and Zhang [24] approximation result is as follows. 

Proof. (Lemma 5.3) Since 4 is a projector from iy 1 (0) onto 1^(0), it follows that on 
each element, 4 is a projector from H 1 ^) onto T\(r), the space of linear polynomials 
over r. Thus, I h p = p, Vp G V\(t), and employing also the stability result in Lemma 5.4 
we have that for < m < k < 2, 

m 

||U - 4u||flm( T ) < ||w-p||fr™( T ) + l|4(p-«)ll# m (T) < C^/l fc_m ||M -p||iffc(s T ), 

fc=0 

where SV is the element support region surrounding r as defined in Lemma 5.4. Employ- 
ing the modified Bramble-Hilbert lemma developed in [10] to estimate the terms of the 
sum gives 

inf \\u - p\\h™(s t ) < Ch k ~ m \u\ H krs T ), 0<m<k<2, 

where due to the assumptions about the domain and the mesh, the constant C depends 
only on the spatial dimension d. Together with the equation above this is 

\\u - I h u\\ H m(r) < Ch k ~ m \u\ H k( ST ) < m < k < 2. 

Since the set 

Q = sup {card{r G 4|r fl5 T / 0}} 
reT h 
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is finite due to the quasi-uniformity and shape-regularity of the mesh, we have finally 
that for < m < k < 2, it holds that 

TETh 

The result for non-integer exponents k and m follows by the usual norm interpolation 
arguments between L 2 {Vt) and H 2 {Vt), which completes the proof. □ 

Lemma 5.3 can be easily extended to the vector case, which provides finally the proof 
of Lemma 3.1. 

Proof. (Lemma 3.1) For u E H^ +Q (fi) = (H^ +a (Q)) d , we have that 

d d 

\\u - I h u\\ 2 L2 = \\ u i - ^Uih^n) < C 2 h 2{1+a) Y klW«(fi)> 
i=i i=i 

where P h denotes the scalar SZ-interpolant applied to Uj. Thus, 



\U - I h u\\ L 2 {n) < Cih 1+a \u\ H i+a( U y 



□ 
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